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Abstract. The drastic decline of the critically endangered northern river terrapin (Batagur baska ) prompted a large-scale 
captive breeding project in Bangladesh and Austria, with the first captive-bred offspring in 2010 . Initially, males and fe¬ 
males were kept together and mated without any system. Flowever, controlled breeding was desired to conserve genetic 
diversity. For revealing relationships among the adult breeding stock and parentages of juveniles, we established a powerful 
genetic marker system using 13 microsatellite loci. Our results indicate that most wild-caught adults of the breeding groups 
are related, suggesting that the wild populations experienced a severe decline long time ago. We develop recommenda¬ 
tions for breeding to preserve a maximum of genetic diversity. In addition, we provide firm genetic evidence for multiple 
paternity and sperm storage in B. baska. Our microsatellite marker system is promising to be useful in breeding projects 
for the other five Batagur species, which are all considered to be Critically Endangered or Endangered. We recommend 
implementing conservation genetic assessments for captive breeding projects of turtles on a broader scale to preserve ge¬ 
netic diversity and to avoid inbreeding. 
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Introduction 

The northern river terrapin Batagur baska (Gray, 1830) is 
a large riverine and estuarine turtle species with a carapace 
length of up to 60 cm (Moll et al. 2009). Like the other 
five Batagur species, the northern river terrapin is charac¬ 
terized by a species-specific, conspicuous breeding colora¬ 
tion in males, which seems to play an important role as 
a prezygotic isolating mechanism (Praschag et al. 2009). 
In B. baska, head and neck of breeding males are black, 
their forelimbs and dorsum of the neck range from light 
orange to rich crimson, and the iris appears yellow (Fig. 1). 
Females are generally distinctly less colourful and are, for 
some species, morphologically difficult to tell apart (An¬ 
derson 1878, Moll et al. 2009, Praschag et al. 2009, un- 
publ. observ.). According to reports of fishermen in the 
Sundarbans, males of B. baska are typically caught in es¬ 
tuaries and along the surrounding coastline, while females 
are found more upstream, especially during the nesting 
season (Moll et al. 2009). 

For a long time, B. baska was thought to be distrib¬ 
uted from northeastern India through the Malay Penin¬ 
sula to Sumatra and Cambodia (Iverson 1992). How¬ 


ever, it was later discovered that the southern and east¬ 
ern populations represent a distinct species, the southern 
river terrapin Batagur affinis (Cantor, 1847) (Praschag 
et al. 2007, 2008, 2009). Batagur baska sensu stricto his¬ 
torically occurred from Odisha and West Bengal (India) 
through Bangladesh to at least the Ayeyarwady and Bago 
river mouths in Myanmar, and possibly to the estuaries of 
the Thanlwin and Sittaung rivers in this country. However, 
B. baska has been extirpated from most of its former range 
with no known wild populations remaining (Praschag et 
al. 2008, Moll et al. 2009). It is now regarded as one of the 
25 most critically endangered turtle species of the world 
(Rhodin et al. 2011) and is considered ecologically extinct 
(Weissenbacher et al. 2015). The species is currently in¬ 
cluded in CITES Appendix I and classified as ‘Critically 
Endangered Aicd’ by the IUCN (Rhodin et al. 2017). 

During the early 1990s, attempts to preserve river terra¬ 
pins now assigned to B. baska sensu stricto failed because 
of political unrest (Moll et al. 2009). Conservation ef¬ 
forts focused on populations in Malaysia, then thought to 
be conspecific. After the recognition of the northern river 
terrapin as a species distinct from B. affinis, efforts were 
resumed, but conservationists needed years to discover a 
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few live terrapins in fish breeding ponds and markets. The 
first individuals were transferred for captive breeding to 
the Turtle Island facility in Graz and the Vienna Zoo, Aus¬ 
tria, where two juveniles hatched in 2010. For large-scale 
breeding, the ‘Project Batagur’ was later launched with 14 
males and six females in the Bhawal National Park (Bang¬ 
ladesh), as a joint initiative of the Bangladesh Forest De¬ 
partment, the Turtle Survival Alliance (TSA), the Vienna 
Zoo, and the IUCN Bangladesh (Weissenbacher et al. 
2015). First breeding occurred there in 2012 and 2013. How¬ 
ever, because of the limited number of breeders, a strategy 
was needed to preserve as much genetic diversity as possi¬ 
ble. To ascertain the relationships among the founders and 
establish such a strategy, we genotyped all northern river 
terrapins from the breeding projects using 13 microsatel¬ 
lite loci. Based on this data, we analysed their relationship 
and developed recommendations for the combination of 
breeding groups. 

Materials and methods 

Sampling 

Currently, 14 males, six females and three wild-caught ju¬ 
veniles are kept in the ‘Project Batagur’ facilities in Bang¬ 


ladesh. Another wild-caught pair and the two first captive- 
bred juveniles (males) live at Vienna Zoo; another male, 
two females and three juveniles are housed at Turtle Is¬ 
land, Graz. Blood and tissue samples of all terrapins were 
studied. An additional male from the project in Bangla¬ 
desh could not be included because it died in 2013, before 
terrapins were sampled. In addition to the adult terrapins, 
23 juveniles from four nests (breeding season 2012) and 61 
juveniles from another four nests (breeding season 2013) 
were sampled (Tables Si and S2). Samples were preserved 
in pure ethanol and kept at -20°C until processing. All ter¬ 
rapins in Bangladesh are microchipped, allowing individ¬ 
ual identification. 


Laboratory procedures 

Total DNA was isolated using the InnuPrep DNA Mini 
Kit for tissues and the InnuPrep Blood DNA Mini Kit for 
blood samples (both kits: Analytik Jena GmbH). Since no 
specific primers for microsatellite loci of any Batagur spe¬ 
cies were available, a wide range of primers developed for 
other turtle species was tested for cross amplification (Sup¬ 
plementary Table S3). The presence of the microsatellites 
in the amplicons was confirmed by sequencing in both di- 



Figure 1. Adult male of Batagur baska in breeding coloration. Photo: Peter Praschag. 
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Table 1. Genetic diversity of studied Batagur baska. n - sample size, nA - alleles per locus, AR - allelic richness, H 0 - observed het¬ 
erozygosity, H e - expected heterozygosity, p - probability, HWE - Hardy-Weinberg equilibrium. 


Locus 

Wild-caught terrapins 

. Private . „ 
n nA .. . AR 

alleles 

H o 

H b 

P 

HWE 

Null 

alleles 

Captive-bred terrapins from Bangladesh 
. Private . „ TT TT 

n nA alleles AR H o H h P 

HWE 

Null 

alleles 

GP19 

33 

8 

2 

7.752 

0.636 

0.772 

0.530 

yes 

no 

83 

6 

- 

5.981 

0.783 

0.712 

0.422 

yes 

no 

GmuB08 

29 

4 

- 

4.000 

0.724 

0.603 

0.643 

yes 

no 

81 

6 

- 

5.877 

0.444 

0.484 

0.743 

yes 

no 

Test21 

33 

5 

1 

4.817 

0.515 

0.622 

0.338 

yes 

no 

83 

4 

- 

4.000 

0.723 

0.637 

0.488 

yes 

no 

TWS190 

32 

3 

1 

2.844 

0.594 

0.496 

0.441 

yes 

no 

83 

2 

- 

2.000 

0.518 

0.398 

0.005 

no 

no 

TWT113 32 

2 

- 

2.000 

0.719 

0.468 

0.002 

no 

no 

78 

2 

- 

2.000 

0.833 

0.489 

0.000 

no 

no 

M01 

32 

3 

1 

2.844 

0.344 

0.446 

0.128 

yes 

no 

80 

2 

- 

2.000 

0.413 

0.418 

1.000 

yes 

no 

M18 

33 

14 

4 

13.262 

0.879 

0.900 

0.440 

yes 

no 

80 

10 

- 

9.981 

0.850 

0.841 

0.000 

no 

no 

M03 

33 

2 

- 

2.000 

0.152 

0.193 

0.298 

yes 

no 

83 

2 

- 

2.000 

0.205 

0.257 

0.078 

yes 

no 

M24 

33 

7 

2 

6.813 

0.727 

0.795 

0.420 

yes 

no 

84 

5 

- 

5.000 

0.714 

0.704 

0.000 

no 

no 

M06 

33 

10 

4 

9.539 

0.788 

0.752 

0.610 

yes 

no 

84 

6 

- 

5.857 

0.702 

0.682 

0.000 

no 

no 

M21 

33 

13 

5 

12.240 

0.879 

0.879 

0.714 

yes 

no 

83 

8 

- 

7.867 

0.867 

0.836 

0.000 

no 

no 

msEo41 

33 

3 

- 

2.969 

0.364 

0.386 

0.777 

yes 

no 

83 

3 

- 

3.000 

0.446 

0.477 

0.024 

no 

no 

GmuD16 27 

14 

5 

14.000 

0.852 

0.891 

0.674 

yes 

no 

72 

9 

- 

9.000 

0.806 

0.843 

0.004 

no 

no 

Mean 

32 

6.769 


6.545 

0.629 

0.631 




81.31 5.000 


4.966 

0.645 

0.602 




SD 


4.423 


4.248 

0.218 

0.214 





2.689 


2.667 

0.196 

0.182 





rections using an ABI 3730 DNA Analyzer (Applied Bio¬ 
systems), and allele size ranges were examined to identify 
informative loci. This resulted in the selection of 13 loci. 

To amplify microsatellite DNA, nine multiplex PCRs 
were performed using a reaction volume of 10 pi contain¬ 
ing 10-20 ng of total DNA, 0.5 units of Taq polymerase 
(Bioron), the buffer recommended by the supplier, 1.5 mM 
MgCl 2 (Bioron), 0.2 mM of each dNTP (Thermo Scientific) 
and 2 pg of bovine serum albumin (Thermo Scientific). In 
addition, each multiplex PCR comprised a specific set of 
primers at a specific concentration. 

Multiplex 1 contained for the loci GmuDsi, GP19 and 
Test56 0.3 pM of each primer; multiplex 2 contained 
GmuBo8 (0.3 pM) and GP55 (1.0 pM); multiplex 3 con¬ 
tained Goag5, GP81 (0.5 pM each) plus GP61 and Test2i 
(0.3 pM each); multiplex4 contained TWL61, TWS190 and 
TWT113 (0.3 pM each); multiplex 5 contained Maucasoi 
(0.3 pM), Maucasi8 (0.5 pM) and Maucas22 (0.6 pM); 
multiplex 6 contained Maucaso3 (0.6 pM) and Maucasi7 
(0.3 pM); multiplex 7 contained Maucaso6 and Maucas24 
(0.3 pM each); multiplex 8 contained Maucasos, Mau- 
casi4 and Maucas2i (0.3 pM each), and multiplex 9 con¬ 
tained GmuDi6 and GmuD55 (0.3 pM each) plus msEo4i 
(0.5 pM). Thermocycling conditions are summarized in 
Supplementary Table S4. Fragment lengths were deter¬ 
mined using an ABI 3730 DNA Analyzer. For this purpose, 
1 pi of the PCR product was diluted in 99 pi water; 1 pi of 
this dilution was mixed with 8.5 pi Hi-Di Formamide (Ap¬ 
plied Biosystems), 0.25 pi water and 0.25 pi GeneScan-600 
LIZ Size Standard (Applied Biosystems). The software 
PEAK SCANNER 1.0 (Life Technologies, Carlsbad, CA) 
was used for scoring fragment lengths. GmuDsi, GP55, 
Maucasos, Maucasi4, Maucas22, Test56, and TWL61 


turned out to be monomorphic, whereas GmuD55, Goag5, 
GP61, GP81, and Maucasi7 had an overrepresentation of 
homozygotes, therefore only the remaining 13 loci were 
used for further analysis. 


Data processing 

Null allele frequencies for the microsatellite loci were ex¬ 
amined using MICRO-CHECKER 2.2.3 (van Ooster- 
hout et al. 2004) and CERVUS 3.0.6 (Kalinowski et al. 
2007). Allelic frequencies, expected and observed hetero¬ 
zygosities were estimated using ARLEQUIN 3.5.1.2 (Ex- 
coffier et al. 2005). The same software was also used to 
test for Hardy Weinberg equilibrium (HWE) and linkage 
equilibrium. Bonferroni corrections were used for both 
calculations. Allelic richness was examined in FSTAT 
2.9.3.2 (Goudet 1995). Statistical significance was exam¬ 
ined using unpaired t tests. 

Maternity, paternity and parent pairs with known sex¬ 
es were determined using the likelihood-based approach 
of CERVUS. Loci with null allele frequencies greater than 
0.05 were excluded to avoid false parental assignments. 
Two levels of confidence were set at 80% (relaxed) and 
95% (strict). Positive LOD scores (i.e. the logarithms of 
the likelihood ratios) were compared to identify the most 
likely parents for each offspring. In 2012 and 2013, five of 
the six females (12654, 12656-12659) and three of the 14 
males (12660,12661,12662) for which samples were avail¬ 
able were kept together for breeding. Thus, using five can¬ 
didate mothers and three candidate fathers, simulations of 
10,000 offspring genotypes were run, each at a sampling 
rate of 100%. The proportion of mistyped loci was set to 
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o.oi. Alternatively, six candidate mothers and 14 candidate 
fathers, with the remaining settings unchanged, were used 
because insemination prior to the separation of the breed¬ 
ing group could have occurred. In addition, relationships 
among terrapins were inferred using ML-RELATE (Ka- 
linowski et al. 2006). In contrast to CERVUS, this Maxi¬ 
mum Likelihood-based software also estimates full sibling 
and half sibling relationships as well as parent-offspring 
(PO) relationships. Moreover, there is a null allele cor¬ 
rection implemented in ML-RELATE, based on a Monte 
Carlo randomization test (Guo & Thompson 1992) and a 
U test (Rousset & Raymond 1995). Thus, estimates with 
and without null allele correction were compared using a 
95% confidence level and 10,000 randomizations. Specif¬ 
ic hypothesis tests were performed using 100,000 simula¬ 
tions. To assess whether a PO relationship is statistically 
more likely than a full sib or half sib relationship, the spe¬ 
cific hypothesis test of ML-RELATE was run. Networks for 
relationships of individual terrapins were visualized using 
PAJEK 3.14 (Batagelj & Mrvar 2003). 


Results 

Even though the captive-bred terrapins showed generally 
less genetic diversity (average allelic richness of 4.97) com¬ 
pared to the wild-caught ones (average allelic richness of 
6.55; Table 1), the differences were statistically not signifi¬ 
cant (t = 1.09, p = 0.29). Also, the number of alleles per lo¬ 
cus was higher in wild-caught terrapins, with their private 
alleles corresponding to those whose alleles were not trans¬ 
mitted to the offspring (Tables 1 and S8). In the genera¬ 
tion, 73.9% of alleles per locus were preserved. Differenc¬ 
es between expected and observed heterozygosities within 
each group were not significant (wild: t = -0.026, p = 0.979; 
captive: t = 0.514, p = 0.612). Furthermore, the data also 
showed no significant difference with respect to observed 
and expected heterozygosities between the groups (H 0 : t = 
-0.117, p = 0.907; H e : t = 0.402, p = 0.691). The deviation 
from HWE in the 1 - generation for some loci (Table 1) is 
expected because the animals originate from a few breed¬ 
ing pairs. 



Figure 2 . Relationship of adult terrapins and hatchlings in Bangladesh, derived from CERVUS and ML-RELATE calculations. Females 
are represented by circles, males by squares, and juveniles by diamonds. Numbers are individual lab codes. Hatchlings from 2012 are 
marked in light grey, hatchlings from 2013 in dark grey. Black lines show PO relations, arrows point from parent towards the offspring. 
Two-headed arrows indicate PO relations between adults with unknown age. For further explanation, see inset. 
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The results of ML-RELATE generally agreed well with 
the assigned parent-offspring relations (PO) revealed by 
CERVUS. Also the assigned full sib (FS) and half sib (HS) 
relationships supported the results of CERVUS. 

Juveniles from four nests hatched both in 2012 and 2013. 
The hatchlings were expected to originate from the five fe¬ 
males (12654,12656-12659) and three males (12660,12661, 
12662) that were then kept together. For each hatchling, 
the parents were identified using CERVUS and ML-RE- 
LATE. In 2012, all hatchlings from the same nest had the 
same mother. Accordingly, female 12657 was the moth¬ 
er of clutch 1 (hatchling 12567) and clutch 2 (hatchlings 
12568, 12570, 12573-12576,12839-12841). For both clutch¬ 
es, a male was revealed as the father (12672) that was not 
kept with the females in 2012. Clutch 3 (hatchlings 12577- 
12587) was laid by female 12659. This clutch was sired by 
two fathers, but only for hatchling 12584 the father could 
be identified (male 12666, again not kept with the females 
in 2012). The father of the remaining hatchlings remains 
unknown. Clutch 4 (hatchlings 12588, 12589) originated 
from female 12654 and from another male (12670) not 
kept with the females in 2012 (Fig. 2; Supplementary Ta¬ 
ble S5). 

For the 2013 season, female 12659 was identified as moth¬ 
er of the 16 hatchlings from clutch 1 (12591,12592, 12594- 
12597, 12600, 12601, 12603, 12605, 12842-12844, 12846- 


12848) and male 12660 as father except for juveniles 12842- 
12844,12847 and 12848, whose father was another uniden¬ 
tified terrapin. Hatchlings from clutch 2 (12845,12.849) re¬ 
sulted from female 12658. The father of hatchling 12845 was 
male 12662, while the father of hatchling 12849 wa s again 
an unsampled male. Female 12657 was the mother of the 20 
hatchlings from clutch 3 (12608-12611,12613-12624,12626, 
12627, 12850, 12851). However, one hatchling (12612) from 
this clutch was assigned to another mother (female 12656), 
which was also revealed as mother of clutch 4 (hatchlings 
12629-12648,12650,12852). In the latter clutch, one juvenile 
(12852) was identified as offspring of female 12657 an d male 
12662. The fathers of the hatchlings from nest 3 are males 
12660, 12662 and 12671, while all hatchlings from nest 4 
most likely originate from male 12660 (Fig. 2; Supplemen¬ 
tary Table S6). As is obvious, many hatchlings have uni¬ 
dentified fathers or fathers (12666,12670-12672) that were 
evidently not kept together with their mothers in 2012 and 
2013. Furthermore, one egg each of clutches 3 and 4 of 2013 
seem to have been inadvertently switched during incuba¬ 
tion (Supplementary Table S6). 

The specific hypothesis test of MF-REFATE suggests 
that most wild-caught terrapins are related in a degree of 
half siblings. Two PO relationships as well as one full sib¬ 
ling relationship were detected among the adult terrapins 
(Fig. 3; Supplementary Table S7). 




Animals Bangladesh 

(used for breeding 2012 and 2013) 

Animals Bangladesh 
(not used for breeding) 

Animals Turtle Island Graz 

Animals Zoo Vienna 

Hatchlings caught in the wild 

Female 


O 


• 


Male 


□ 


■ 


Juvenile 




♦ 

O 


—*— PO relationship 

- FS relationship 

-HS relationship 


Figure 3. Relationship of wild-caught terrapins as well as of animals bred in Vienna Zoo and Turtle Island, derived from CERVUS and 
ML-RELATE calculations. Numbers are individual lab codes. The red line shows full sib relations; dashed lines, half sib, grandparent- 
grandchild, uncle/aunt-nice/nephew, or first cousin relations. For further explanation, see Figure 2 and inset. 
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Table 2. Recommended mating combinations for captive adult Batagur baska (parent-offspring and half sib relationships excluded; the 
only full sib relationship found refers to two males). * - combinations that produced offspring in 2012 and 2013. 



12656 

12657 

12658 

12659 

12854 (Vienna) 
12858 (Graz) 
12862 (Graz) 


n<r 

no* 

no 

no* 

no 

no 

yes 


yes 

yes 

yes 

yes 

yes 

yes 

yes 


yes 

no* 

no* 

yes 

no 

yes 

yes 


yes 

no 

yes 

yes 

no 

yes 

no 


yes 

yes 

yes 

no 

yes 

yes 

no 


yes 

yes 

yes 

yes 

yes 

yes 

yes 


yes 

yes 

yes 

no 

yes 

yes 

yes 


yes 

yes 

yes 

no 

yes 

yes 

no 


yes 

yes 

yes 

yes 

yes 

yes 

yes 


yes 

yes 

yes 

yes 

yes 

yes 

no 


yes 

yes 

yes 

yes 

yes 

yes 

yes 


yes 

yes 

yes 

yes 

yes 

yes 

yes 


yes 

no* 

yes 

yes 

yes 

yes 

yes 


yes 

yes 

yes 

yes 

yes 

yes 

yes 


no 

yes 

yes 

yes 

yes 

yes 

yes 


yes 

yes 

yes 

yes 

yes 

yes 

yes 


Discussion and conclusions 

Our study is one of the relatively few examples for the ap¬ 
plication of conservation genetics in a captive breeding 
programme for reptiles (see the review in Witzenberger 
& Hochkirch 2011). These authors surveyed 188 studies 
using conservation genetic approaches for breeding pro¬ 
grammes. The majority of studies applied microsatellites, 
with a clear focus on mammals and birds, and only 9% 
were referring to reptiles. However, turtles belong to the 
most endangered of any of the major groups of vertebrates 
and their status is paralleled only by the similarly endan¬ 
gered primates (Rhodin et al. 2017). Among turtles, the 
northern river terrapin (Batagur baska ) is one of the most 
critically endangered species of the world (Rhodin et al. 
2011). For the conservation of this ecologically extinct spe¬ 
cies, and many other turtle species, captive breeding plays 
a key role (Williams & Osentoski 2007, Weissenba- 
cher et al. 2015). Yet, most captive breeding projects are 
not monitored genetically and only recently attention has 
been paid to the pedigree or relatedness of breeders using 
conservation genetic approaches ((^ilingir et al. 2017 for 
the closely related B. trivittata ), even though Williams & 
Osentoski (2007) recommended the implementation of 
such measures one decade ago. 

The breeding projects for B. baska in Bangladesh and 
Austria aim to preserve as much as possible the current ge¬ 
netic diversity. To achieve this goal, we developed a power¬ 
ful micro satellite marker system to assess individual rela¬ 
tionships. This marker system allows estimating the rela¬ 
tionship of wild-caught and captive-bred individuals, and 
thus constructing pedigrees to plan future breeding. 

In Bangladesh, 84 juveniles hatched during the seasons 
2012 and 2013. We studied these terrapins, the wild-caught 
breeders, and all terrapins from the Austrian breeding col¬ 
onies. Our results are promising. The first captive-bred 
generation retained approximately 75% of the genetic di¬ 


versity of the wild-caught founder individuals (Table 1), 
even though no breeding strategy had been implement¬ 
ed yet. We are confident that with introducing a breeding 
strategy combining less related or unrelated individuals, 
the percentage of preserved genetic diversity can be sig¬ 
nificantly increased. 

According to our results, most wild-caught terrapins 
are related (Fig. 3). Given the long life expectancy of tur¬ 
tles, this situation suggests that the wild population ex¬ 
perienced a severe decline long ago. The few survivors 
are largely related at the level of half sibs (or first cousins, 
aunts/uncles-nices/nephews or grandparents-grandchil- 
dren). For conservation purposes, the present genetic di¬ 
versity should be preserved to the greatest extent possible. 
Thus, the reproduction of closely related terrapins has to 
be avoided. Based on our data, we suggest that selective 
breeding should be implemented, and some mating com¬ 
binations be avoided (Table 2). Moreover, we recommend 
that in following seasons the males 12661, 12663-12665, 
12667-12669, and 12673, which did not contribute to previ¬ 
ous offspring, should be used for breeding. In addition, the 
adults from the two breeding colonies in Austria should 
also be involved for breeding in Bangladesh. The highest 
priorities have those unrelated individuals listed in Table 2 
that have not bred before. 

A remarkable finding for the colony in Bangladesh was 
that some captive males contributed to offspring that must 
have mated with the females prior to the reproductive sea¬ 
sons of the present study. Moreover, we could not identi¬ 
fy the fathers of some juveniles (Fig. 2). Even if the latter 
should be all offspring from the same unsampled male that 
died in 2013 in Bangladesh, it was not kept together with 
the females in 2012 and 2013. This provides unambiguous 
evidence for sperm storage in B. baska. Also females of 
other turtle species are known to store sperm in their ovi¬ 
ducts, so that clutches can be fertilized long after the last 
mating (Gist & Jones 1989, Pearse & Avise 2001, Pearse 
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et al. 2001, Roques et al. 2004, 2006, Johnston & Rand 
2006). Our results also provide evidence for multiple pa¬ 
ternity in B. baska (Fig. 2; Tables S5 and S6), a phenom¬ 
enon supported by sperm storage and thought to increase 
fertility rate (Olsson et al. 1996) and genetic diversity (An- 
dersson 1994). 

Our microsatellite marker system developed for B. bas¬ 
ka also promises to be informative for the five other spe¬ 
cies of the genus Batagur. All species are considered by the 
last IUCN Red List assessment as either Critically Endan¬ 
gered ( B. affinis, as part of B. baska; B. borneensis; B. ka- 
chuga; B. trivittata) or Endangered ( B. dhongoka ) (Rhodin 
et al. 2017), and for many taxa captive breeding projects 
are underway. Microsatellite primers usually perform well 
within, and often beyond, the same turtle genus (unpubl. 
observ.; see also Schwarz et al. 2003, King & Julian 2004, 
Vamberger et al. 2011, 2017, Tiedemann et al. 2014). Thus, 
we expect to have provided a most useful tool for preserv¬ 
ing genetic diversity in breeding projects for Batagur spe¬ 
cies in general. 

Compared to the ddRAD-Seq approach to generate 
SNPs used by Qilingir et al. (2017) for B. trivittata, our 
microsatellites are an easy and affordable method to gen¬ 
erate an informative marker set that does not need high- 
quality DNA. Our approach can be easily repeated by any 
basic molecular genetic laboratory. 

Finally, we explicitly recommend implementing con¬ 
servation genetic assessments for other captive breeding 
projects of turtles to preserve a maximum of genetic diver¬ 
sity and to avoid inbreeding depression. 
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